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ABSTRACT 


The main objective of the research presented in 
this thesis is to calibrate Indo Swedish Traffic Simulation 
Model for free speed situation. The road submodel has 
been modified accordingly. Extensive validation of the free 
flow simulation have been carried out. 

In the first phase major effort has been devoted in the 

study of free speed distributions of different types of 

vehicle. Free speed model, as a submodel of the main traffic 

simulation model, considers the effect of roadway parameters 

(i.e. width, curvature and roughness) on the basic desired 

speed (BDS). Models have been developed to compute the free 

speed of vehicles as function of roadwidth, curvature and 

roughness. The median basic desired speed, is reduced to 

Om 

a new median speed V , where V, is a function of road 

Im Im 

width. Then V is reduced to a new median speed , which 
Im 2m 

is a function of both of the median speed without the 

horizontal curve, V , and the radius of curvature. Now V_ 

Im 2m 

is reduced to a new median speed which is a function of 

3m 

surface roughness and median speed V„ . After the median 

2m 

speed V has been computed, a resulting new speed 
3m 



distribution is calculated. This is obtained by transforming 
the basic desired speed distribution with the median 
value V to the new median speed V. and at the same 
time rotating about V so that the dispersion of the 

w 

distribution decreases. A transformation measure, the Q 
value, indicates how far the basic desired distribution must 
be rotated about the new median value V, , is computed for 
each vehicle type. 

The road submodel has been modified to 
calculate the median speed and transformation measure, Q 
value, for each vehicle type. The free flow simulation model 
is than validated for each vehicle type for different road 
geometry in plain and rolling terrain for the single lane, 
intermediate lane and two lane roads with varying horizontal 
curvature and roughness. The predicted free speed are found 
to be very close to those obtained in the field, thus proving 
the validity of the calibration of free flow simulation 
model. The model thus developed is capable of predicting 
accurately the free speed distribution of vehicles on any 
road stretches thus filling one of the gap existed in the 


Indo-Swedish Traffic Simulation Model. 



CHAPTER I 


INTRODUCTION 


1 . 1 General 

Transportation contributes to the economic, industrial, 
social and cultural development of any country. Among all 
modes of transportation, road transport is the most easily 
accessible mode available to people. The road network could 
alone serve the remotest villages of the country. Road 
transport has recorded a phenomenal growth in the wake of all 
round national development in the country. Hence there is a 
need for enhancing and upgrading the road network. This 
causes substantial pressure on the planners. In order to 
select the most beneficial scheme from among a number of 
alternatives, it is necessary to have an economic appraisal 
of various alternatives. Rational decisions will be required 
before undertaking such improvements. This necessitates the 
availability of reliable information regarding free speed, 
time delay, fuel consumption, accident rates and cost of 


maintenance and spare parts. Given the roadway and its 



geometric features, the above mentioned factors form the 


function of speed and quality of flow on these highways. 

1.2 Need for Road Development Planning 

In the above framework there is an urgent need to 
develop rational tools to evaluate alternative investment 
programs. The investments on the proposed schemes are of very 
high magnitude. Budget and resource constraints forces us to 
closely examine the various costs and benefits associated 
with the investment programs such as construction costs, 
maintenance costs, benefits in the terms of user costs. Such 
a framework would enable the planner to select a project 
judiciously among the various investment programs under 
consideration. 

1.3 Indian Road and Traffic Condi tioiis 


Indian 

roads 

comprise primarily 

of 

single lane. 

intermediate 

lane. 

two 

lane and 

four 

lane divided 

carriageways. 

The 

road 

with pavement 

width 3.5m are 


designated as single lane roads, while the roads of 5.5m wide 
are generally known as intermediate lane roads. 7m wide roads 



are the two lane highways. The four lane divided carriageway 


roads are 14m to 15m wide, provided with median barrier. 
In general all the roads are provided with earthen shoulders 
on either side. Some times part of the shoulders of the 
single lane road has been improved with brick paving to 
increase the effective width of the roads to facilitate the 
ease of movement on these roads. The brick paving vary in 
width from about 0.5m to 1.6m. However, their condition is 
generally poor due to usage by heavy vehicular traffic. Since 
the strength of the shoulders is much lower, they deteriorate 
much faster when subjected to wheel heavierloading. Some 
times the difference in levels is as much as 10cm, at the 
interface between the pavement and the brick shoulders. 

The road traffic on Indian roads is highly 
heterogeneous, being constituted by vehicle types as fast as 
Maruti cars and as slow as bullock carts. The fast moving 
vehicles include light motor vehicles, trucks and buses, 
trailer combinations and several categories of two wheelers. 
Amongst the slow moving vehicles, bicycles, cycle rickshaws 
and animal drawn vehicles dominate the scene, 
It has been observed that operating speeds of these vehicles 


are quite low even at low and medium traffic conditions. This 



may be attributed due to poor surface conditions, road 
geometric, over loading of vehicles, vehicle conditions etc. 
The traffic on some sections of national highways has grown 
to such an extent that proposals for widening the two lane 
highways to four lanes are under consideration. The 
Government policy to develop some locations as industrial 
areas and the acute traffic congestion on many links has 
created the need for the construction of roads with superior 
facilities such as expressways. In addition to this, there 
are numerous proposals to widen some of the existing single 
lane roads to intermediate and two lane highways. 

1.4 Earlier Work Done in India 

Realizing the need for setting up certain design 
standards, evaluation of various project schemes and the 
necessity to create a database on various road conditions, 
traffic characteristics, a project Road User Cost Study 
CRUCS) was sponsored by the Government of India and the 
International Bank for Reconstruction and Development (World 
Bank). In order to evaluate the alternative investment 
policies, the road user cost study was carried out. The road 


user cost study in India dealt with only this component of 



the total transportation cost. 


Evaluation of traffic flow characteristics can be made 
in two ways: through analytical studies and simulation 
modeling. The simulation modeling project was taken up by the 
Indian Institute of Technology, Kanpur(IITK) in 
collaboration with the Swedish National Road and Traffic 
Research Institute( VTI) . The study developed simulation 
models for single lane, intermediate lane, two and four lane 
conditions with extensive calibration and validation. The 
simulation models developed for the Indian traffic conditions 
for single lane, intermediate lane, and two lane roads 
consider the crossing of vehicles in opposite directions. 
These models were developed with Swedish model as their 
basis. 

1.5 Indo Swedish Traffic Simulation model 

This model is the modified version of the VTI model, 
originally developed in Sweden, to suit the road and traffic 
conditions prevailing in the Indian environment. On the basis 
of its application the model has been divided into two parts. 


One part simulates the speed of the freely moving vehicle 



along the road "free flow traffic model" and another part 
simulates the interaction between the individual vehicles in 
the vehicle streams "interaction model". 

1.5.1 Free Flow Traffic Model 

The road stretch is represented in the model as a series of 
geometrically homogeneous road blocks. Within each block the 
factors of road width, radius of curvature, speed limits and 
gradients are kept constant. Each driver in the traffic 
simulation model is allotted a speed from the basic desired 
speed distribution(BDS) . A driver unimpeded by traffic 
interactions chooses to travel at maximum speed on a wide, 
straight, horizontal road with no speed restriction and this 
is termed as basic desired speedCBDS) 

The free block speed adopted by a vehiclein a 
particular road block is dependent on the median speed of 
the road block, transformation coeff icient(Q value) and BDS 
of the vehicle. Apart from the block speed, the slope, if 
any, is also assigned to the block. Each vehicle has also 
been allotted a power mass(P/m) ratio, as well as air and 


rolling resistance, which together decides the ability to 



reach or maintain this block speed on the particular slope. 


If the vehicle has a speed lesser than the block speed, when 
entering a block, or it cannot maintain its block speed due 
to the slope, its speed is calculated using Newton’s equation 
for momentum. 

The resulting speed is thus determined from a desired 
ability situation resulting in a speed profile along the road 
in which corrections have been made to the block speed with 
regard to P/m ratio of the vehicle to allow for the vertical 
profile of the road and the air and rolling resistances. 

1.5.2 Interaction Model 

The free vehicle model calculates how free vehicle 
adopts its speed to road width, alignment, gradient and speed 
limits. The interaction model takes into account the way in 
which individual vehicles are influenced by surrounding 
traffic. The speed of the vehicle in the traffic stream is 
influenced by the volume, composition and speed distribution 
of the traffic. 


In addition to the speed and P/m ratio, each vehicle 



is allocated a time and road coordinate for starting 


a 


direction of travel and an initial entry speed. This permits 
the faster vehicles to catchup with slower vehicles along the 
road and leads to overtaking or following. The overtaking has 
been divided into flying, accelerating and passing. If 
overtaking’ starts as soon as one vehicle catches up with 
another vehicle, flying overtaking take place. If overtaking 
occurs after following, an accelerated overtaking is carried 
out. If the road has a wide shoulder of good quality 
overtaking can also take place by passing, which means that 
the leading vehicle moves with certain probability on to the 
shoulders. When following occurs, the vehicle follows the 
preceding vehicle with a constant time headway depending on 
the preceding vehicle type. 

In addition to the free movement, vehicle meet each 
other in opposing directions, decelerate to safe crossing 
speeds and cross each other. Two opposing vehicles are 
interacting with each other due to the width of the road on 
single lane and intermediate lane roads. On these roads 
vehicles have to yield space and slow down, if the two 
vehicles involved in crossing are wide. Under Indian 
conditions it is more important to consider crossing. 
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since speeds are calculated 


based on roadwidth, shoulder 


type and its quality and vehicles type involved in 
crossing. 

1.5.3 Pi'ogram System for Traffic Simulation 

The simulation program comprises the heart of a 
program system as shown in the Figure 1.1. The simulation 
model is a model system permits computation of the following 
detailed events: 

a) Simulation of the behaviour of free flow traffic in 
both directions, i.e., vehicle speed 
prof ile’sdependence on road attributes, speed limits 
in force on vehicle attributes, but excluding the 
interdependence of other vehicles in the traffic 
stream. 

b) Simulation of the behaviour of vehicles which also takes 
into account of the interdependence of surrounding vehicles 
in their own and opposite direction. 


1 . 5. 4 Road submodel 



The objective of this model is to represent the road as 


a series of geometrically homogeneous road blocks and to 

create a speed profile along the road for each free moving 

vehicle. Thus a road block will be of constant road 

width, curve radius, gradient, speed limit and roughness. 

The speed profile depends on road geometry and speed limit. 

The median speed for each block V , is a function of the 

3m 

distribution of BDS, road width, horizontal curve, and the 

roughness. In addition a transformation measure, the Q 

value for each road block is calculated, which 

indicates how far the BDS distribution must be rotated about 

the median speed V to obtain the free speed 

3m 

distribution. 

1.5.5. Traffic Submodel 

The traffic submodel describes the traffic 
characteristics provided as data input to the simulation 
model. Traffic in the simulation model is represented by 
number of driver vehicle(DV) units. A DV unit contains the 
characteristics to model the driver behaviour as well as the 


vehicle parameters. The traffic is made up of fast as well as 



slow 

moving 

vehicles. 

Vehicles vary widely 

in 

their 

size. 

power 

to mass 

ratio, operating speed 

and 

other 


characteristics such as vehicle age and condition. The 
attributes of a DV unit are 

a) Identification number. 

b) Vehicle type. 

c) P/m ratio 

d) Basic desired speed. 

e) Origin. 

f) Destination. 

g) Entry time 

h) Entry speed. 


Every DV unit in the simulation model is uniquely 
represented by an index called identification number. 
Associated with each DV unit are the vehicle parameters, 
vehicle type and P/m ratioCp-value) . The p-value is a 
significant factor contributing to the performance capability 
of the vehicles. Lower p-value limits the ability of a 
vehicle to climb the grades thereby affecting the capacities 


of the road stretches considerably. 



1.6 Problem Statement and Object Of Study 


It has been observed that different vehicle types have 
different BDS distributions and the variation between the BDS 
distributions is appreciable. Specifically in the case of 
Indian traffic, it is necessary to consider the median block 
speed and Q value separately for each vehicle type. Thus it 
has been decided to calibrate the Q value for ten different 
vehicle types under different road condition and geometry. 
Once the calibration of Q is done, the necessary modification 
in the road submodel is carried out and an extensive model 
validation for free speed simulation is performed. 

Therefore, the main objective is to calibrate Q value 
for ten different vehicle types from the field data collected 
by the Central Road Research InstituteCCRRI) Delhi. The 
second step is the modification of the road submodel and 
validation of the simulation model for free speed condition. 

1.7 Organization of the Thesis 


The thesis is presented in six chapters. Chapter II 



presents the detailed description of the Indo Swedish Road 
Traffic Simulation model, traffic submodel and road submodel. 
A brief description of the models are present in this 
chapter. 

Chapter III explains in detail the model for free speed 
prediction. Calibration of various parameters . involving VTI 
free vehicle model are explained in this chapter. This 
chapter also explains the derivation of polynomials for 

free speed prediction. 

Chapter IV explains the free speed modeling and 
calibration of Q values. Free speed distributions generated 
using the calibrated Q values are also shown in this chapter. 

Chapter V explains the various inputs to the road 
submodel, experiments with the road submodel and validation 
of the model. The comparison between the field and simulated 
data are done in this chapter. 

Chapter VI presents the summary of the present work. 


scope of the study and conclusion drawn from the study. 



CHAPTER II 


REVIEW OF THE LITERATURE 


2. 1 Gener al 

The detailed description of the Indo Swedish Traffic 
Simulation Model for narrow, intermediate and two lane roads 
is presented in Section 2.2. The traffic behaviour in the 
simulation model considering overtakings are given in 
Section 2.2.3. The detailed description of road submodel 
covering from calculation of median speed to the 
transformation coef f icient(Q value) are given in section 2.3 
The generation of traffic from calibrated distributions using 
traffic submodel is presented in Section 2.4. 

2.2 Indo Swedish Traffic Simulation Model 

A stochastic discrete event simulation model was 
developed by the VTI for the two lane bi-directional traffic 
system during the period 1965 to 1977. It has a long history 
of calibration and validation over a number of road 
stretches in Sweden. Subsequently, the model has been used 


for traffic analysis in providing auxiliary lanes in Finland 



and United Kingdom. While RUCS study was initiated in India 


the need to develop a simulation model for the Indian road 
and traffic conditions was felt. Accordingly, an evaluation 
of several models existing before 1980 led to the conclusion 
that the best available and reliable model was the one 
developed by VTI. This model formed the basis for work done 
subsequently at the Indian Institute of Technology, Kanpur 
in collaboration with the scientists of VTI. The version of 
the model is currently known as the Indo-Swedish Traffic 
Simulation Model. 

The traffic and road condition prevailing in India 
are too complex to model by simple approaches. It is 
necessary to incorporate the heterogeneity of traffic and the 
hosts of road widths along with shoulders condition. The 
basic structure of the VTI model is such that it had built 
in features that allowed the restructuring of the model for 
Indian conditions. The program and the data structures are 
based on the Jackson Structured Programming( JSP) technique 
and has been programmed in one of the highest level of 
languages - SIMULA67 The main program principle of SIMULA is 
pseudo parallel execution which can be regarded as a further 
development of event control of dynamic sequences. The model 


has been calibrated and validated extensively for the single 



2.2.1 General Program [>esign 


The program mainly consists of two processes, which also 
contain the data and procedures. They are 

1. Process class generator process 

2. Process class vehicle 

The vehicle generator process creates the driver 
vehicle objects and allots the individual drivei — vehicle 
attributes. Here they are also allotted their traffic 
attributes. Parameters which define the vehicle are 
identification number, basic desired speed and power mass 
ratio. Parameters giving their traffic attributes are origin 
of the vehicle, destination of the vehicle, entry time and 
entry speed. The vehicle generator process also activates 
vehicles at their starting times. 

The vehicle process describes all possibilities for action 
that a particular vehicle can have, for example “drive as 
freely moving vehicle" , "follow another vehicle”, "overtake 
the vehicle in front" , "change to another track" etc. The 
freedom of choice covers mainly the set of alternatives that 


any vehicle can have, and the actions are assumed to occur 



momentarily at calculated times. At each event the model data 


is updated and a particular event generated from among the 
possible event types. A note of the predicted event is 
inserted inserted chronologically and logically in a two way 
linkedlist and the events are then executed in this order. 

The ordinary cycle for an arbitary vehicle is: 

1. Predict the time of the next event - PREDICTNEXTEVTIME 

2. Wait for the predicted time - HOLD 

3. Move the vehicle in time and space - DRIVE 

During the phase one of the cycle 

PREDICTNEXTEVTIMECpredict next event time), 

PREDBLOCKBORDERTIME(predict time of passage of next 
blockborder) , AVERSPC speed from preceding event to next 
event) and PREDICTBLBORDERSP(predict speed at passage of next 
block border) are calculated in their corresponding 
procedures. During phase three of the cycle the procedure 
DRIVE updates the attributes LOCALTIMEC time of preceding 
event), LOCALCOORD( road coordinate at preceding event), 
AVERSP, PREDBLBORDERTIME(predicted time of passage of next 
road block border), and PREDECTBLBORDERSP(predicted speed at 


passage of next road block border). 



During the phase two of the cycle it may happen 
that an adjacent vehicle interacts with the current vehicle 
with the result the predicted time of the next event for this 
vehicle is shown incorrect, since it will occur earlier. Its 
ordinary cycle consequently isinterrupted and the current 
vehicle considers that a surprise has occurred through the 
SURPRISE procedure. A prediction of new earlier event 
time for the vehicle is then made. 

In the program the stretch of the road consists of a 
sequence of consecutive road block objects and a sight 
distance function in each direction of travel. Each road 
block object is homogeneous with regard to the roadwidth, 
slope, horizontal curvature, roughness, speed limit and 
overtaking restriction. The road block is represented in the 
program as an object from the Link class Roadblock. The 
object has the following important attributes : 

-RCO coordinate of the roadblock 

-RBLLENGTH road block length 
-ROADWIDTH road width class 

-LANE occurrence of auxiliary lane/ lateral space 

-RI slope 

-RCURVE radius of curvature 

-ROUGHNESS roughness of the road 


-RQ 


Q value with respect to normal speed 



-IR 


block median speed 


RQ RQ 

-DVQ IR - (median speed of the vehicle) 

~CR block cross median speed 

~CQ Q value Mith respect to crossing speed 

CQ CQ 

-CVQ CR - (median crossing speed of the vehicle) 

CQ, CVQ and CR are considered for accounting the 
traffic behaviour on narrow roads. 

2.2.2 Representation of Traffic Behavior 

Figure 2.1 shows the representation of lanes across the 
road width. A vehicle unimpeded by other traffic normally 
chooses to travel in lane 2. Thus lane 2 is used when driving 
under normal conditions whereas lane 3 is used while yielding 
to a faster vehicle from behind, and lane 1 is used while 
overtaking a slow vehicle. When vehicle P is caught up behind 


by a 

faster 

vehicle, Q, in lane 

2 then P would 

move on 

to 

lane 

3 with 

certain probability 

so that, the 

catching 

up 


vehicle could pass P in lane 2 (Fig 2.2). 

When the catching up vehicle Q is not allowed to 
pass in its normal lane, i.e., when a paved shoulder is not 
present, or when P rejects to move on to the shoulder(if 
shoulder exists), Q would look for a gap in lane 1 and starts 
overtaking without having to reduce its speed. This is termed 



as "flying overtaking". Figure 2.3 shows the steps involved 
in flying overtaking. When the gap available in lane 1 is not 
acceptable by the caught up vehicle Q, then it prepares to 
follow P at a safe distance. When an opportunity to overtake 
is accepted then Q moves to lane 1 for an accelerated 
overtaking. Thus, lane 1 is normally used for overtaking 
maneuvers, lane 2 when unconstrained by other traffic or 
while following a vehicle in front, and lane 3 to allow a 
faster vehicle catching up from behind to pass in lane 2. 
This logic with some modifications has been applied for 
single lane, intermediate lane and two lane bi-directional 
movement simulation. 

Each vehicle has been allotted a "head" and a "tail". 
The driving rule in this model tries to maintain a minimum 
gap between any two vehicles, in the same lane. In order to 
ensure that this gap is maintained, the head length and tail 
length are defined for each vehicle as follows. 


Headlength 


(V - V )' 

forward follow 


2 * Deacc 


If V, .<V. _ 

forward follow 

= 0 otherwise 


( 2 . 1 ) 


Taillength = Th^ ,, * V, , 

follow follow 

where V, . is the speed of the preceding vehicle, 

forward 


V,, , , is the speed of the following vehicle, 
follow 



Th, ,, is the time headway of the following 
follow 

vehicle. 

Deacc is the retardation rate. 

Head length therefore indicates the distance required 
to slow down with the given retardation rate to the speed of 
the preceding vehicle. The tail length is chosen so that 
vehicles in queue are separated by given time intervals. 

As shown in Figure 2.1 a roadst retch is divided into 
homogeneous blocks. Each block has median speed IR, Q value 
and RQ in addition to other attributes. In the case of narrow 
roads median cross velocity(CR) , Q value for cross velocity 
and CVQ for each block are also present. The DVQ and CVQ are 
calculated as explained in section 2.2.1. 

When a vehicle enters a block its free block 
speed and in the case of narrow roads free block and its free 
cross speeds are calculated by using the following formulae: 

Free block speed = (VON® + (2.2) 

CQ 1 / CQ 

Free cross speed = (CON + CVQ) (2.3) 

where VON is the basic desired speed of the vehicle 

CON is the median crossing speed of the vehicle 



A numerical integration of the power equation over time is 
applied to determine the speed and time of reaching the block 
border. The procedure continues till the free block speed has 
been exceeded or block border has been reached (Fig. 2.4) 

2.2.3 Probabilities of Overtaking 

The drivers desire to perform an overtaking operation with 
a given free sight distance is determined in the model by a 
set of stochastic functions. The probability of accepting a 
free sight distance of a given length when opportunity for 
overtaking is offered, depends on different overtaking 
situations as described below: 

Type of overtaking 1. flying 

2. accelerating 

Sight limitation 1. oncoming vehicle 

2 natural obstruction 

Overtaken vehicle 1. vehicle type 

2. speed of the overtaken vehicle 
Roadwidth 1. hard shoulder > 2m, corresponding 

to a paved road width of 11m. 

2. others 


The probabilities are given by 



p(x) = 0 


if x<= Sim. 


= a(x-Sl)/(S2-Sl) if Sl< X < S2m (2.4) 

= others if S2 <= x 

where x is a free sight distance and a,Sl,S2 are 
calibration constants. 

In addition, if the vehicle is in platoon, the 
probability is reduced as follows : 

P . = p*kCplace in platoon - 1) (2.5) 

red 

where P ^ is the reduced probability 
red 

p is obtained from the equation 2.4 
k is model constant at present it is set to 0.75 

2, 3 Road Submodel 

This model creates the road as a series of homogeneous 
road blocks and the variation of sight distance along the 
road, as input to the simulation program. In the simulation 
model each vehicle has an essential characteristic, i.e. its 
BDS. This is the speed at which the driver would chose to 
travel on a wide straight horizontal road with no speed 
limit, i,e, "free speed". Brodin and Carlsson (1983)CRef .3) 
described free speed as a function of road geometry and BDS, 


which formed the basis for the VTI model and Indo-Swedish 



Traffic Simulation Models . The road submodel described by 
Brodin and Carlsson is explained in this chapter. This is 
modified to a more general form and is explained in the next 
Chapter. 

It is assumed that for a straight horizontal road having 
width of 12m. or more, there is aBDS distribution and from 
that vehicle chooses to travel at a speed known as its BDS. 
The vehicle is prevented from maintaining its BDS by the 
following factors: 

* a road width less than 12m. 

* curves with a radius of less than 1000m. 

* speed limit 

* upward gradients. 

Effect of Roadwidth: 

For a roadwidth under 12m. the median basic desired 
speed is reduced to a median speed assumed 

that roads with over 7m paved are built with a carriage way 
7m wide, the remaining width consists of two hard shoulders. 
For a roadwidth less than 7m it is assumed that no hard 


shoulders are provided. 



The median speed V is given by: 

im 


1 

V 


Im 



a 

<v 


a 

475 


( 2 . 6 ) 


where v. is the road width in m 
b 

a is the calibration constant 

is the median speed for 7m. road width. 
Om 


Effect of Horizontal Curves: 


If a road includes a horizontal curve, the median speed 


V, must be reduced taking into account the horizontal 
Im 


radius. Curves with mean radius r > 1000m. do not effect the 


speed V . If r <= 1000m. the median speed V_ on the curve is 
Im 2m 


obtained by the following expression, where speeds V„ and 

2m 


V are given in m/s. 
Im 


2m 


SQRT(C-~)^+ b (-- 0.001)) 

V r 

Im 


(2.7) 


where r is the mean radius in meter 


b is the calibration constant. 


V is the median speed on straight road. 
Im 


From the equation it is observed that, entry speed V 



small radius curves constant bcon has the gratest effect. 


Effect of surface roughness: 


For surface roughness, the median speed V, is obtained 

3m 

with the following expression: 


V_ = V :»: ( ) (2.8) 

3m 2m 


where Rg is the surface roughness 

c is the calibration constants 


After the median speed, V, has been calculated, the 

3m 


resulting distribution V must be determined. A 

w 


transformation measure, Q, indicates as to how far the basic 


desired speed distribution must be rotated about the median 


speed V, . The Q value is a function of the median speeds 
3m 


, V„ and V_ , The free speed distribution V^, using 
Om Im 2m 3m 3 


the Q value transformation is determined as follows: 




V 


Q 

3m 


(2.9) 



in basic desired speed and free speed distributions 
respectively. 

If Q = 1 then the above equation results in a purely 
parallel shift indicating a constant reduction in speed for 
fast movipg as well slow moving vehicles. 

However, when Q <1 the free speed distribution, V , is 
rotated in anticlockwise direction showing that a driver with 
a higher BDS reduces his speed more than a driver with a 
lower BDS when influenced by speed reduction factor. It must 
be noted that the smaller the value of Q, the larger will be 
the rotation. 

Each of the factors of the road width, horizontal curve 

and speed limit has its own isolated measure of rotation, 

termed q. 

This factors are: 

q^ for road width 

q for horizontal curve 
2 

q for speed limit 
3 ^ 

All values of q. < 1 

1 

The total measure of transformation, Q, is a function of q. 

with median speed V ,V„ and V, included as weighing factors 

Im 2m 3m 



Q 


ql kl + q2 k2 + q3 k3 
kl + k2 4 k3 


( 2 . 10 ) 


where k = V - V 

1 Om Im 

'<3 = ''sJ 


Hence Q is the weighted mean of q ,q and q 

X ^ w 


Speed profile for the median vehicle: 

The road to be simulated is described as a series of 
homogeneous blocks. Within each block, the factors of road 
width, curve radius, surface roughness and also gradient are 
constant. Thus V is constant within each block. If any road 

O 

factor is altered, a discontinuity in the speed profile 
results. Consequently, retardations and accelerations between 
these blocks are treated as follows: 

If the median speed V, in block k is greater than 

3m 

the median speed V, in block k4i, the vehicle is retarded in 
block k until the next median speed is reached at the block 
limit. A new block limit is created where the retardation 


begins. 



block k 



block k+1 


V 


3 


On the other hand, if speed V, is less in block k 

3in 

than k+1, the vehicle accelerates in block k+1 to the true 

median speed V, . This acceleration takes place with a 

3m 

special speed model, the p-model, described in the next 
section. 



Fig. 2. 6 

At each retardation, but not acceleration, a new block 
is created, termed the influence block. The length of this 
influence block is obtained as follows: 


K,k+1 

are the 

block number, 

3m, 

. and 
k 


. .are 
,k+l 

the 

corresponding 

desired 

median speeds, calculated 

in m/s. 

In 

the influence 

block. 

the retardation 

r is 

constant. 

The 

length L of the influence block is then 

given 

by 

the motion 



r 


m/s 


( 2 . 11 ) 


V 


2 

3m , k4 1 


2L 


V 


2 

3m, k 


2 

The r value is taken as 0.5 m/s . This 
corresponds approximately to the retardation a vehicle 
receives when the driver releases the accelerator pedal. This 
gives 


L 


3m, k 


V 


2 

3m, k+1 


( 2 . 12 ) 


According to the magnitude of L, the length of the 
influence block is chosen as follows: 


'-inf 


L (L >= 25m) 

25 (13 <= L < 25m) 
0 (L < 13m) 


(2.13) 


If L < 13m, a retardation is obtained within 
block k+1, but this error is insignificant. 

k 


2.4 Traffic Subn^xlel 


This model assigns the vehicle characteristics(identity 
number, basic desired speed and P/m ratio) as well as traffic 
characteristics(s tar ting point, starting time and 
destination). This model was originally developed by Clive 
Gilliams of U.K Department of Transport(Ref .2) . The current 



model is the modified version. The distribution of vehicles 
among the various classes (i,e,, cars, trucks .scooters etc.) 
is to be specified in the input data. The procedure for 
generating the vehicle as well as traffic parameters is 
explained below. 

Identification Number: This is trivially generated by 
incrementing each time a vehicle is generated to enter the 
simulated road stretch at the specified coordinate. 

Vehicle type: A distribution of vehicle types is specified 
for the anticipated traffic composition. The vehicle type of 
an individual vehicle is obtained from the distribution. 

Basic desired speed: The vehicles are uniformly distributed 
over the basic desired speed distribution for each vehicle 
type as specified in the input. It generates the 8DS from the 
given distribution. 

Power to mass ratio: The power equation can be obtained as 
follows. 

Consider a vehicle on the upgrade. Figure 2.6 shows the 
various forces acting on the vehicle. The force equation 
gives 

= f - F^- -rag* sin(i) 


m a 



C2.14) 


where m = mass of the vehicle in Kg. 

2 

a = acceleration in m/s 
F = tractive force in Newtons 


Fj^= air resistance in Newtons 

F = rolling resistance in Newtons 
r 

g = gravitational acceleration in ms 


-2 


i = gradient in m/km 


V = vehicle speed in m/s 


t = time in s 


The tractive force at the wheel can be written as 


F = P/V 


(2.15) 


where P is the power available, measured at the wheels 
in Watts and V is the speed in ms 
The air resistance is expressed as 




(2.16) 


where is an air resistance coefficient in kg/m' 

2 

A is the frontal area in m 


The rolling resistance F^ is expressed as 

F = m cos(i)* (C , + C _V) 
r rl r2 


(2.17) 



c 

rl 


= tn CC , 4 C 
rl 

and C are 
r2 


r2''> 

rolling 


resistance coefficients. 


The gravitational force is expressed as 
m g sin(i) ~ m g i 


(2.18) 


Substituting the above expressions in force equation, 
we get 


dv P 2 

m --- = -- - C AV - m(C 4C ) - mgi 
dt V L crl cr2 


dividing by m we get 


dv _ p 
dt ' ~V 



CCri+ -gi 


(2.19) 


P 

where p = is the power/weight ratio used 

m 

for the vehicle (p-value) 


The equation is solved numerically by conversion to a 
difference equation 


dv .. Ay_ _ p 
dt At “ V 


C A 

- (C 4 C V) - gi (2.20) 
m rl r2 


Speed is 


V 


n4i 


thus calculated as follows 

V 4 A V 
n n 


n = 


0 , 1 , 2 . 



2 


AV = At[ --- 
n ^ V 
n 



^ ‘'rl* 


gi] (2.21) 


And distance covered is 


S = S + A S (2.22) 

n+1 n n ^ 

As = 0.5 . At C V 4 V , ) 

n n n4 1 


Time stage At is chosen as 


At 


0.015 



(2.23) 


before iteration is started. 

Where is the entry speed in the block 

a^ is the acceleration at this speed 
0 


a - -B !r-_ V ' 

O V m 0 

0 


■- (C 4 C ^ V^) - gi 
rl r2 0 


(2.24) 


So that the time stage At is not too large for small a^, it 
is limited to the same value as for travelling half the block 
length 


At 


block length 


(2.25) 


max 2 V 

In this way, at least two iterations are obtained in each 


block. 



Entry time: The entry time for each vehicle entering 
the simulated road stretch is generated using the 

Shchul ’s(1955) composite time headway model. The usual form 

of the composite distribution is 

fCt) = (1- oi)* g(t) + a * h(t) (2.26) 

where f(t) is the probability density function of the 

composite headway distribution, g(t) is the probability 

density function of the headway distribution for free moving 

vehicle, h(t) is the probability density function for 

constrained vehicles, and a is the proportion of the 

constrained vehicles. Branston (1976) who analyzed the 

constrained time headway on the two lane roads in Indiana, 

USA, found that a lognormal distribution fitted the data with 

a mean of 2 seconds and a standard deviation of logarithms of 

constrained headways o of 0.4. The time headway distribution 

c 

of free moving vehicles is determined by platoon length, n. 
The total headway for a platoon length n is the sum of the 
one headway for the free moving platoon leader and (n-1) mean 
headways for the constrained vehicles. The mean time headway 
for the stream can be expressed as 

% 

+ (P - 1) M (2.27) 

q T c 



where (J is the mean of the platoon size, q is the 
traffic flow in veh/h and /u and jj are the mean time headways 

T c 

for free and constrained vehicles respectively. A two 
parameter distribution due to Miller(1974) has been used for 
the generation of platoon size. The distribution of a 
platoon, size r is defined as 


pr 


(s + r - 1) 

Cm + s + r + 1) 


if r > 1 


(2.28) 


(m + 1) 

(m + s + 2) 


if r=l 


(2.29) 


Miller (1965) found that jU may be estimated by using the 
equation 


Ai(est) = 0.58 + 1.58 * z 


(2.30) 


where z = 


0.1 * q 
X (1 - q /j^) 


(2.31.) 


where X is the overtaking rate for constrained 
vehicles. Now, the problem of estimating the M^ is reduced to 
that of estimating X for a given road traffic flow. Miller 
(1963) used data collected on a single straight two lane road 


in Sweden to produce the estimate 



the oncoming flow is assumed to be the traffic flow given as 
input, q. 


Entry speed : For a free vehicle the entry speed is taken 

to be equal to 0.85*BDS . In the case of a constrained 

vehicles, the minimum BDS among the platoon is calculated 

(BDS . ) and 0.85 * BDS . , , is assigned as the entry 

minpl minpl 

speed of all the vehicles in the platoon. However, all these 
speeds are constantly reassessed once the vehicle enters the 
road stretch to be simulated. Hence entry speeds assigned are 
comparatively less important. 

Detailed description of the existing traffic simulation 

model, road submodel and traffic submodel is presented in 
this chapter. The models for free speed prediction are 
presented in the next chapter. Chapter IV explains the 
calibration of Q values for different vehicle types. Chapter 
V explains the various experiments carried out for validation 
of the simulation model. 
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(a) Vehicle Q catching up with P 
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PIQ. 2.2 Pssing manoeuvre described in Indo-Swedish 








(a) Vehicle Q catching up with P 



FIG. 2.3 


Flying overtaking procedure described in 
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CHAPTER III 


MODEL FOR FREE SPEED PREDICTION 


3. 1 Gener al 

Free Speed model, as a sub model of the main Traffic Simulation 
model, considers the effect of roadway parameters (i,e, width, 
curvature, roughness, gradient etc ) on the basic desired 
speed(BDS). Thus study of BDS and Free Speed of different types 
of vehicles becomes essential, as a prime input for the traffic 
simulation model. 

A driver may not always be able to drive his vehicle at 
the basic desired speed (BDS) because of certain constraints like 
inadequate road geometry, riding quality of pavement and traffic 
flow conditions. The speed of the operation of any vehicle is 
governed by factors like roadway, traffic and other factors such 
as driver behaviour and vehicle characteristics (make and model). 

The spot speed of any vehicle at any particular 
location depends upon a number of factors such as ; 

* Road factors (road geometry, road surface 
characteristics etc). 

* Traffic factors (Traffic intensity, traffic 
composition etc). 



TABLE NO.3.1 FRAMEWORK FOR FREE SPEED MEASUREMENT 
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* Vehicle factors (make and model, maintenance 


condition etc). 

* Human factors ( driver behaviour etc ). 


3.2 Field Data Presentation 

The requirement of the free speed study includes the effect 
of road width, road roughness and road curvature in plain terrain 
and rolling terrain. Taking all these parameters into 
consideration the framework that spells out the cells 
representing the road stretches with the specific road type, road 
roughness and radius of curvature where the free speed have been 
recorded for each of the vehicle types. 

In order to assign a particular number to each of these cells 
a system of numbering was evolved by CRRI . Each of the site in 
a cell has been numbered as FSCxyz 

FSC represents Free Speed on Curved sections. 

X represents road width 

1 for single lane 

2 for intermediate lane 

3 for two lane 

y represents road roughness 

1 for Rg less than 3000 mm/Km 


2 for Rg between 3000 and 6000 mm/Km 



3 for Rg more than 6000 mm/Km 


z represents radius of curvature and terrain 

1 and 6 for Rc less than 150m in plain terrain 
and rolling terrain respectively 

2 and 7 for Rc between 150m and 300m in plain 
and rolling terrain respectively 

3 and 8 for Rc between 300m and 450m in plain 
and rolling terrain respectively 

4 and 9 for Rc more than 450m in plain terrain 
and rolling terrain respectively 

5 and 0 for straight approach section in plain 
terrain and rolling terrain respectively. 


3.^. Planning of the Study 

The speed of only those vehicles are included which 
are being driven under free flow conditions and not 
influenced by either the road sight distance or because of the 
traffic conditions. Road stretches satisfying the requirements of 
each and every individual cell of the matrix framework as shown 
in the Table no. 3.1 were identified. 

Free speed of different types of vehicles were recorded with 
the help of radar speedmeters suitably located on the straight 
and curved sections of the selected road stretch with appropriate 



road geometry. 


Road geometry and the road roughness were measured 
for each individual site with the help of instruments 
like gyroscope, gradometer and car mounted bump integrator 
units etc by the Central Road Research Institute (CRRI). 

The speed distribution thus obtained for various vehicle 
types in each cell of the above referred matrix was analyzed to 
test for the normal distribution through chi-square test as 
reported by the CRRI. 

The median speed and the standard deviations as 
calculated for the observed set of data for each road type in 
plain or rolling terrain, horizontal curvature and road 
roughness are tabulated as per CRRI identifiers. 

3.4 Free Vehicle Model of Indo-Swedish Traffic Simulation Model 

The aim of the free vehicle model is to create a speed 
profile along the road for each free vehicle. This profile 
depends upon road geometry. The calculation is done in several 
steps according to the approach given below. 

3.4.1 Roadwldth Model 

According to this model, if the road width is <= 7m, 

median speed V is obtained according to the following 



are given in m/s. 


expression, where speeds V and V 

Om Im 


1 


V 

Im 


1 


V 

Om 


+ 


a 

o;— ^xs") 


a 

2.5 


(3.1) 


of 


where is the BDS in m/s 

Om 

y is the road width in m 
b 

a is a calibration constant 
for all vehicle types. 

Fig. 3.4.1 shows the median speed 
roadwidth for Ambassador car. 


with 



a value of 0.0042 


as a function 


3.4.2 Model for Horizontal curve 


If the road includes a horizontal curve, the median speed 

V must be reduced taking into account the horizontal radius. 
Im 

Curves with a mean radius r>= 1000m do not affect speed V . If 

Im 

V <= 1000m the median speed V on the curve is obtained with 

Im 2m 

the following expression, where speeds V and V are given in 

2m Im 

m/s. 



1 

SQRT((-^ )^ + b (C-^ ) - 0.001)) 
Im 


(3.2) 


where r is the mean radius in meters 
b is a calibration constant 

Fig. 3. 4. 2(a) to Fig. 3.4.2(c) shows median speed as a 

49 



function of curve radius r. 


3. 4* 3 Sur face Roughness Model 


For a given roughness Rg, the median speed V, is obtained 

3m 

with the following expression where V and V are speeds in m/s 

3m 2m 


V 


3m 



( 


-c*Rg 

e 


) 


(3.3) 


where c is the calibration constant 

Rg is the surface roughness in mm/km. 

Fig.3.4.3Ca) to Fig.3.4.3Cd) shows the median speed V 

3m 

as a function of free speed V_ and surface roughness Rg. 

2m 


Since the free vehicle model of Indo-Swedish Traffic 
Simulation Model considers only four vehicle types and has the 
same basic desired speed distribution for all the vehicle types, 
the number of parameters to be calibrated for free speed 
prediction are vary few. But considering the heterogeneity of the 
Indian road traffic, it is desired to assign separate basic 
desired speed distribution for each vehicle type. Thus the number 
of parameters that need to be calibrated is quite large. Moreover 
it is observed that, free vehicle model slightly over estimates 
the free speed when applies to Indian traffic conditions. Hence 
it has been decided to derive a polynomial fit for free speed 


prediction instead of free vehicle model, which gives more 
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accurate prediction. The various steps involved in polynomial fit 
for free speed prediction has been explained in the subsequent 
sections below. 

3. S Poiynoniials for Free Speed Prediction 


In order to predict the free speed for each and every cell in 
the above mentioned matrix framework, which is latter used in 
the main simulation program, number of polynomials are formulated 
using the observed data. The polynomials formulated for different 
road geometry are as follows: 

1. Polynomial fit for roadwidth CV^). (order =2) 

2. Polynomial fit for curvature (V^): 

This has three different polynomials 

(a) For single lane (SL) road . (0rder=2) 

(b) For intermediate lane road (IL) (order=3) 

(C) For two lane road (TL) (0rder=3) 

3. Polynomial fit for roughness (V^): 

This has following different types of polynomial: 

(a) For single lane road: 

i. For curvature <= 300m (order=2) 

ii. For 300m <curvature<= 450m (order=2) 

iii. For 450m <curvature<=1000m (order=2) 


iv. For straight road Corder=2) 



(b) For intermediate lane road: 

i. For curvature <=150m Corder=2) 

ii. For 150m <curvatu re<=300m Corder=2) 

iii. For 300m <curvature<=450m Corder=2) 

iv. For 450m <curvature<=1000m (order=2) 
V. For straight road Corder=2) 

(c) For two lane road: 

i. For curvature<=150m Corder=2) 

ii. for 150m <curvature<=300m (order=2) 

iii. For 300m <curvature<=450m (order=2) 

iv. For 450m <curvaturB<= 1000m (order=2) 
V. For straight road(order=2) 


3.5.1 Polynomial for Roadwidth 


The form of the polynomial that have been fitted here is as 
follows 

2 

Y= a +a. .X + a„.X 
0 1 2 

where Y= Free speed of a particular vehicle. 

X= Road width in meters. 

a . a . a are coefficients of the polynomial. 

0 1 2 

Since the order of the polynomial is two, three sets of 
equations are required to compute the value of the 


f 



computation of coefficients are 



a 

+ 


-X 

+ 

a . 

X 

1 

0 


1 

1 


2 




+ 


-X 

+ 

a . 

X 

2 

0 


1 

2 


2 


Y = 

a 

+ 

a 

.X 


a . 

X 

3 

0 


1 

3 

+ 

3 



2 
2 
2 

Where » Y2 ’''^3 free speed of a vehicle on single lane, 

intermediate lane and two lane roads respectively with minimum 
surface roughness on a tangent alignment (i,e straight road). 
The corresponding cell numbers are FSS 115, FSS 215, and FSS315. 
^2*^2’ ^3 “ Roadwidth corresponding to single lane, 

intermediate lane and two lane roads respectively (i.e.3.5m, 5.5m 
and 7.0m). The computed coefficients for different vehicles are 
shown in the Table 3.5.1, 


3.5.2 Polynomial for Radius of Curvature 

In order to predict the free speed of a vehicle on a curved 

segment of a road of radius of curvature X meter for a 

particular road width a polynomial of the form 

2 

Y = a + a .X + a_.X for single lane road 

0 1 2 

2 3 

Y = a^ + a^.X + a^.X + for two and three 

lane road, 

were formulated. Here Y indicates free speed of a vehicle on a 
curved road of radius of curvature X meter. In order to compute 


the coefficients of the polynomial three sets of equations for 



single lane and four sets of equations for intermediate and two 
lane roads are solved. The set of equations that are solved 
for single lane are: 

2 

2 

> 

2 


^0 * ^2-Xl 


^3= *1-X3 ^ ^2-X3 


Where ''^j[ speeds corresponding to the cell 



number 

FSC 

112, FSC 

113, 

FSC 

114 respectively. 

X ,X„ X, are 

1 2. 3 

the radius 

of curvature 

150m, 

375m, 

and 725m 

respectively. 

Since 

observations 

have 

been 

made 

in three 


different groups of horizontal curve (less than 300m, between 
300m and 450m, and 450m and 100m) the mid values of each group 
has been taken. 


For intermediate lane and two lane roads the sets of equations 
that are solved for computing the coefficients are 

Xr ^0 ^-Xl " ^2-Xl^ X 

2 3 

Y = a + a .X„ + a„.X„ + a^.X^ 

20 12 22 32 

2 3 

Y,= a^ + a .X, + a^.X, + a_.X_ 

30 13 23 33 

2 3 

Y = a^ + a, .X ■ + a -X + a .X 

401424 34 

Where Y ,Y ,Y ,Y are speeds corresponding to the cell number 
12 3 4 

FSC 211, FSC 212, FSC 213, FSC 214 respectively for intermediate 
lane and FSC 311, FSC 312, FSC 313, FSC 314 respectively for two 
lane road. The values of X^, X^, X^ and X^ taken are 75m, 225m, 


375m and 725m respectively( i.e. mid values 


of less than 


150m, 150m to 300m, 300m to 450m, and 450m to 1000m). The 



3*5*3 F’ol ynonii All for Surface Roughness 


In order to predict the free speed of a vehicle for any given 

roughness, curvature and roadwidth, a polynomial of the form 

Y= a + a .X + a .X^ 

0 1 2 

is formulated. Here Y represents free speed of a vehicle for a 
given roadwidth, radius of curvature and of roughness X mm/Km. 

coefficients of the polynomial. In order to 
compute the coefficients of the polynomial, three sets of 
equations for a given roadwidth and radius of curvature are 
solved. The equations that are solved are as follows 

Yi= a^.X^H. a^.x/ 

^2= V *1-V 

Y_= a.+ a .X_+ a .X ^ 

3 0 1 3 3 3 

Where Y ,Y ,Y are the speed corresponding to the cell number FSS 
12 3 

115, FSS 125, FSS 135 for single lane straight road. And X ,X ,X 

JL jC ■C) 

are the values of surface roughness 1500 mm/Km, 4500 mm/Km and 
7500 mm/Km respectively(i .e. mid values of less than 3000 mm/km, 
between 3000 mm/km and 6000 mm/km, and more than 6000 mm/km). The 
computed values of the coefficients are shown in the Table 3.5.3. 

The predictions have been made with the polynomials as 
derived above, and the predicted values were found to be very 
close to the actual values. The polynomials were subsequently 
implemented in the road submodel, and the results obtained by the 



explains the different steps involved in calibration of Q Value 
for different vehicle types. 



TABLE 3.5.1 COEFFICIENTS OF THE POLYNOMIAL FOR ROAD WIDTH EFFECT 


VEH. TYPE 

^0 


®2 

AMBASSADOR 

56.24 

0.94 

-9.50X10~^ 

FIAT 

59.12 

1.02 

-1.42X10~^ 

MARUTI 

63.99 

1 .28 

-1 .43X10~^ 

JEEP 

51.77 

0.84 

-9.53X10~'^ 

SCOOTER 

45.95 

0.72 

-1.90X10~^ 

MATADOR 

50.46 

0.69 

-9.53X10~'* 

BUS 

58.98 

0.76 

9.51X10"“^ 

LCV 

50.95 

0.66 

1.42X10~^ 

TRUCK 

54.72 

0.82 

-1.42X10~^ 

TRACTOR 

25.87 

0.32 

-4.76X10'^^ 



TABLE 3.5.3Ca) COEFFICIENTS OF THE POLYNOMIAL FOR ROUGHNESS EFFECT 
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TABLE 3.5.3Cc) COEFFICIENTS OF THE POLYNOMIAL FOR ROUGHNESS EFFECT 
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CHAPTER IV 


FREE SPEED MODELING 


4. 1 Gener al 


Free speed (FS) of a vehicle is the speed it can 
actually maintain on a road stretch in the absence of any 
traffic interaction. In case of ideal road condition FS 
equals BDS (V^)- However, if road geometry, including surface 
characteristics imposes any restrictions, FS would be lower 
than the BDS, the extent of reduction depending upon the type 
and magnitude of the constraints. 

The main road factors influencing travel speed, 
resulting in FS value lower than BDS are identified as: 

1 . Road width 

2. Horizontal curvature 

3. Surface roughness 


4. 2 Model for Transformation of Basic Desired Speed Distribution 

We know that basic desired speed is reduced to 

median speed V when influenced by road width, horizontal 
3rri 

curvature and surface roughness. 

It is now necessary to create a distribution of the 

speed about the median value V^, termed asthe desired speed 



termed as Q value. 


4.2.1 Importance of Q-Value in Siniulation Program 

In the main simulation program for each road 
block and for each vehicle free block speed is calculated 
using the following equation. 

Q 1/Q 

freeblock speed = (VON + DVQ) 

Where VON is the basic desired speed. This basic 

desired speed for a vehicle is obtained from the 
BDS distribution of the vehicle from the traffic 
file, 

Q Q 

DVQ = V - VOM 
3m 

VOM = Median basic desired speed distribution. 

V = Median block speed which is obtained after 
3m 

considering the roadwidth, horizontal 
curvature and pavement surface roughness. 

Q is the measure of rotation which shows how far the free 
speed distribution is shifted because of the roadwidth, 
curvature and pavement surface roughness. 



4.2.2 Importance of Q-Value in Traffic Generation 


The same Q value is used for determining the basic 
desired speed of the vehicle in traffic generation 
programme for the field traffic data. The procedure used for 
calculating the BDS for the field traffic is as follows. 

The basic desired speed of the vehicle from the given BDS 
distribution is given by 

BDS = v(j,i)® + DVQ(j,i)^^® 

Where vCj,i) is the basic desired speed of the i th type 
vehicle at j th percentage. Percentage 'j* is drawn 
randomly from the inbuilt procedure. 


Q is the rotation factor for i th vehicle. 

DVQCj.i) is calculated by using the equation 

DVQ(j.i) = - V3M ° 

where V is the median basic desired speed and V3M 
Om 

is the observed median speed at the site. 

Hence the value of Q is an important factor in 
determining the free block speed for each vehicle in the 


simulation program and 


in determining 


the BDS from 



the spot speed in traffic generation 


program from the 


field data. 

4.2.3 Del er mi nat i on of Q- Val ue 


The Q value is a function of the median speeds , V, , 

Om Im 


V , and . Using the free speed distribution V_, the Q 
2m 3m 3 

value is determined as follows. 



Where V and V . are speeds at an arbitrary percentile in 
Oi 3i 

basic desired speed and free speed distributions 
respectively . 


If Q = 1 then the above equation results in a purely parallel 
shift indicating a constant reduction in speed for fast moving as 


well as slow moving vehicles as shown in the Fig. 4. 2.1 


However, when Q < 1 thefree speed distribution, V^, is 
rotated in anticlockwise direction showing that a driver with 
a higher BDS reduces his speed more than a driver with a 
lower BDS when influenced by speed reducing factors as 
shown in the Fig. 4. 2. 2 . It must be noted that the smaller 
the value of Q, the larger will be the rotation. Each roadway 
factor namely the roadwidth, horizontal curvature and 
pavement surface roughness has its own isolated measure of 




Distribution 



Fig.4-Z«l Effect of Road Geometry on BDS Distribution 

Uniform reduction at high as well as low speed 



4.2*4 Determination of q. Values; 


A calibration constant q < 1 such that V, is 

i - i 


approximately 


distributed as V minus a shift, d. , 

o t 


which is dependent only on 


the factor i. Thus 


describes the effect of the factor i on 

q. q. 

i.e. v/ = V V d. 

l o 


the shape of V . 

o 


Hence for every sample v^ from V^, v^ from 


"•i "*1 

V. = v^ + d 
X 0 


q.-l qj-1 

Differentiating: q. v. dv. = q. v^ dv^ 


1-q. dv 

[ — 9— ] ^ = [ ] 

■■ V ■* dv. ■* 

i 1 


1 - q. 


log [ 


log [ 


s constant 


From the above equation q^ value can be obtained. In the 
case of VTI model these constants were calibrated for the 
Swedish road and traffic conditions and they are: 


0.5 for road width. 



= -0.8 for horizontal curve. 

= -2.0 for speed limit. 

These values were used for all the vehicle typesCi.e. four 
types of vehicle) and they had adopted only one BDS 
for all vehicular traffic on Swedish roads. 

■ In the VTI model these constants were calibrated 
without differentiating the roadwidth as they have mostly 7m 
or wider roads. But in the case of Indian roads roadwidth 
vary from 3.5m wide single lane roads to 7.0 m wide two 

lane roads.lt is interesting to observe the rotations 

with respect to road width for each vehicle type. Hence 

the rotation factors for each vehicle type, on each type 

of roadwidth for the effect of road width, horizontal 

curve and surface roughness, q, can be obtained. 
Where i = 1,2,3 (i-e. effect of roadwidth, radius of curve, 
surface roughness.), j = 1 to number of vehicle 

types. The approximate procedure for obtaining the q. . 

is as follows. 

1. Determine the BDS distribution with 2% increment in 
the cumulative distribution for each vehicle type. 

2. Determine the free speed distribution at 3.5m, 5.5m, 

and 7.0m wide roads with the effect of roadwidth, 

I 

horizontal curve, and roughness can be obtained with 2% 



1 nc r©rn©rit in "t-h© cufnul 3 tiv© 
vehicl© typ©. 


distribution 


for 


each 


3. The value of q at roadwidth w and for vehicle type 
j can be calculated as follows. 


1 - 


q 


log(dv . /dv .) 
qj ij 


log(v . / V ,) 


4. Similarly these distributions are to be calibrated 

at a particular horizontal curve on each of the road 

width class, (i.e. at 3.5m, 5.5m and 7.0m). From these 

distributions as mentioned above q and q are to be 

2 j 3 j 

calculated. 


The following example shows the determination of q value 

'•J 

for given distributions of V .and V .. Fig. 4. 2. 3 shows the 

OJ 

assumed distribution of V , and V .. Two tangents are drawn 

ij oj 

at 50 percentile on both the distributions as shown in 


the Fig. . The values of dv . ,dv v . and v . are 

Oj IJ oj Ij 

measured as shown in the Figure. The value of q . is 

U 

calculated using the above equation. Distributions of 


li i 

V ' . V " were calculated using q value, 

ij oj iJ 

Fig. 4. 2. 4 shows the variation of these distributions. 
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distributl ons 


Table 4.1 Values of V 


of oj 




and V 


Hj 


IJ 


Cumulative V , 

Percentage 

V 


o 

> 

OJ 1 J 

5 

68.0 

76.0 

0 . 000907 

0.000754 

-0.00015 

10 

72.5 

81.0 

0.000816 

0.000679 

-0.00013 

15 

75.5 

85.0 

0.000763 

0.000626 

-0.00013 

20 

78.0 

88.5 

0 . 000722 

0.000586 

-0.00013 

25 

80.5 

92.0 

0 . 000686 

0.000549 

-0.00013 

30 

82.5 

94.5 

0 . 000658 

0.000525 

-0.00013 

35 

84.0 

97.0 

0.000639 

0.000503 

-0.00013 

40 

85.5 

99.5 

0.000620 

0.000482 

-0.00013 

45 

87.5 

101.5 

0.000597 

0.000466 

-0.00013 

50 

88.5 

103.5 

0 . 000586 

0.000452 

-0.00013 

55 

90.5 

105.5 

0.000564 

0.000437 

-0.00012 

60 

91.5 

107.5 

0.000554 

0.000424 

-0.00013 

65 

92.5 

109.5 

0 . 000544 

0.000411 

-0.00013 

70 

94.0 

112.0 

0.000530 

0.000396 

-0.00013 

75 

96.0 

115.0 

0.000512 

0.000379 

-0.00013 

80 

98.0 

117.5 

0.000494 

0 . 000366 

-0.00012 

85 

100.0 

121.0 

0.000478 

0.000348 

-0.00012 

90 

102.5 

125.0 

0.000459 

0 . 000330 

-0.00012 

95 

107.5 

132.5 

0 . 000424 

0.000300 

-0.00012 

100 

115.0 

147.5 

0.000379 

0.000251 

-0.00012 


equally 


spaced. Table 4.1 shows the values used for 


oj 


distributions 

and 

V . '^-V . From 

these 



Oj ij 


it can be concluded 

that 

the value 

of q . 




IJ 


describes parallel shift between V 


oj 


V . distributions, 
ij 


4 . 2.5 


Calculation of Values: 


The total measure 
q (j = 1, 2, 3 and 
the median speed 
factors . 


of transformation Q. is a function of 


i = 1 


1, 

to number of vehicle types) with 

V V .included as weighing 

2m,i 3m,i 


Q = 

LV 


q, . k . + q_ . k . + q 
1,1 1,1 2,1 2,1 3,1 

k - . + k „ . + k _ . 

1,1 2,1 3,1 


1,1 1,1 Om, 1 im, 1 


k„ . = a .( V, -- .) 

2,1 2,1 lm,i 2m , 1 


k, . = « .CV^ - - V, .) 
3,1 3,i2m,i 3m , 1 


The value of Q. 

t 

q. 

V '■ approximately 
3i 


describes the 
distributed as 


shape of the V . such that 

31. 

Ol. 

L 

V . minus shift d, which 
o i 


is 


dependent on all these factors. Hence the weighing 



Hence we propose the following q values are to calibrated. 

i-j 

For narrow lane roads and for each vehicle type: 

1. q for road width effect for different road width range. 

2. q for radius of curve effect for different radius of curve. 

3. q for roughness effectt for different roughness limits. 

3j 

The above mentioned procedure have been 
tried out for the calibration of values, but the results 
obtained were not satisfactory. Since the rotation parameter 
are calculated on the basis of BDS data(mean and standard 
deviation) and median block speed, the difference in means 
gives a measure of lateral shift and the difference in 
standard deviation gives the measure of rotation. The more 
the individual SD values, more is the dispersion and less the 
rotation. Again, each median block speed SD must be less than 
the SD of its BDS distribution because more the vehicle speed 
more is the SD- The data supplied by the CRRI were not 
consistent with respect to the above mentioned factors. The 
inconsistency in the data may be because of the following 

reasons: 

i) The data that have been collected were not from 
and horizontal curvature and roughness were 


the same location 



ii) The vehicles observed on different location 

were not the same vehicle. 

Thus it is decided to calibrate the values by 

trial and error i.e. we started with an assumed q. value and 

1 

a distribution is predicted. The predicted distribution is 
then compared with the actual distribution and their sum of 
squares of the difference is computed. The q. value for which 
sum of the squares of the difference is minimum is adopted. 
The subsequent section explains the various steps involved in 
calibration. 

4. 3 Free speed study 

Free speed model, as a sub-model of the main Traffic 
Simulation model, considers the effect of roadway parameters 
(i.e. roadwidth, curvature, roughness etc.) on the basic 
desired speed(BDS). Thus study of BDS and free speed of 
different types of vehicle becomes essential, as a prime 
input for the traffic simulation model. 

As per the accepted definition in Traffic 
Engineering, basic desired speed, at which a driver would 
desire to drive his vehicle on a 2 - lane rural highway with 
wide hard shoulders, provided there are no obstructions 
due to traffic and the speed is not restricted by geometric 

t 
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cell no FSS 315 satisfy all the conditions stated above and 
hence taken as basic desired speed. 

4.3.1 Objectives 

To study the distribution of the basic desired speed 
(V^) for each vehicle type. 

To estimate the median speed (V ) and rotation 

O 

parameter(Q) for each vehicle type 

4. 3. 2 Scope of study 

The scope of the study is confined to cover the 
following parameters 

Type of roads : 

The following type of roads are covered 

- Single lane road 

- Intermediate lane road 

- Two lane road 

Type of terrains: 

The study will cover following two terrains 

- Plain terrain 

- Rolling terrain 

Plain terrain for the purpose of study are 
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meters pr Km. and rolling terrain for total rise and fall 
between 15 to 40 meters per Km length of the road. 

Type of Vehicles: 

Vehicles to be covered are classified into following 
sub groups 

- Light motor vehicles 

* Ambassador car 

* Premier Padmini 

* Maruti cars and vans 

* All makes of jeeps 

* Two wheelers (scooters and motorcycles) 

- Heavy motor vehicles 

* Matadors and mini buses 
^ All makes of buses 

* Light commercial vehicles 

* Trucks up to 12 tonne capacity 

- Slow moving vehicle 

* Agricultural tractors 

4.4 Calibration of q^-Values 

AS it is stated earlier that, each roadway factor 
namely roadwidth, horizontal curve and pavement surface 
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roadwidth effect, i=2 for curvature effect, i=3 
for roughness effect). Thus it is necessary to calibrate 
vales in order to estimate the total measure of rotation Q. 
The subsequent section explains the the various steps 
involved in calibration of q. values. 


4 . 4.1 


Calibration of 


value 


ql represents rotation factor for road width alone. This is 

calibrated by using the transformation 

V -V 

Oi 3i Om 3m 


or 


V . = W ^-(V 
3i 0 Om 3m 


1/ql 


Where - Median value of 8DS 
Om 

V = Median value of free speed 
3m 

V . =BDS at an arbitrary percentile 
Oi 

V .=Free speed at an arbitrary percentile 
3i 

For calibrating ql values for different vehicle types 
the mean and standard deviation given in the cells FSS315, 
FSS215 and FSS115 as shown in the Table 4.4.1 are taken and a 
two percentile distribution is generated. Using the 
transformation as given above values are predicted for the 
cells FSS115 and FSS215. The predicted values are then 
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of the difference for both the cells are computed. Initially 
ql value is assumed to be 1.0 and it is then decremented 
by 0.02 after each iteration. The value of ql for which sum 
of the squares of the difference is minimum is adopted as ql. 
Fig. 4. 4.1 shows the plot between ql and total sum of squares 
of the difference. The procedure is repeated for each 
vehicle types. 

4.4.2 Calibration of value 


q2 represents the rotation factor for horizontal curvature 
alone. This is also calibrated by using the transformation 

V . = [V 

Zx Ox Om 3m 

For calibrating q2 values for different vehicle types the 
mean and standard deviation given in the cells FSC311, 
FSC312,FSC313,FSC314 as shown in the Table4.4.1 are taken and 
two percentile distribution is generated for each of 
thecel Is. Using the transformation stated above V^. values are 
predicted for the cells FSC311,FSC312, FSC313, FSC314. The 

predicted values are then compared with the field data and 
their total sum of the square of the difference for all 
the cells are computed. Here also initially q2 value is 
assumed to be 1.0 and it is then decremented by 0.02 after 

each iteration. The value of the q2 tor which total sue of 

adopted. Fig. 4. 4. 2 shows the plot 


the squares is minimum is 
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procedure is repeated for all 


vehicle types. 


4.4.3 Calibration of value 


q3 represents rotation factor for surface roughness 
alone. This is also done by using the transformation 

V =tv 

3i ^ Oi ^ Om 3m 

For calibrating q3 values for each of the vehicle types 

the mean and standard deviation given in the cells FSS315, 

FSS325 and FSS335 are taken as shown in the Table 4.4.1 and 

two percentile distribution is generated. Using the 

transformation as given above, V values are predicted for 

oX 

the cells FSS325 and FSS335 - The predicted values are then 
compared with the field data and their sum of squares of 
the differences are computed. Here also an initial value of 
q3 is assumed as 1.0 and it is then decremented by 0.02 
after each iteration. The q3 value for which sum of the 

I 

squares of the difference is minimum is adopted. Fig. 4. 4. 3 
shows the plot between q3 and sum of squares of the 

difference. The procedure is repeated for each vehicle types. 


4.5 Calibration of a value 

The total measure of rotation of Q is a function of ql, 
q2,and q3 with median speeds V^^included as 
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Where = MCV^^- V^^) 

^ = “(''le- ''zJ 

^3 = 

V , V, , V^ and V_ are the median speeds. 
Om Im 2m 3m 


4.5.1 Value of a. 


Here assumed value of ai is taken as 1 for all types of 
vehicles. 


4. 5. a Calibration of 

When only road width and curvature effects are taken into 
consideration the expression for Q reduces to 

''iJ ‘’1 * ''2m) ''2 

In the above expression only a2 and Q values are unknown. 
To calibrate the ca values for different vehicle types we 
assumed a positive value for «2 and Q value is computed from 
the above expression. Once Q is calculated V^. values for the 



cells FSC112, FSC113, FSC114, FSC211, FSC212, FSC213. FSC214 
FSC311 , FSC312, FSC313, and FSC314 as shown in the Table 

4.4.2 are predicted. The predicted values are compared with 
the field values and their total sum of squares of the 
differences is computed. The Qt2 value is incremented by 0.02 
after each iteration. A graph is plotted between ot2 values vs 
their corresponding total sum of the squares of the 
differences as shown in the Fig. 4. 5. 2. From the plot it is 
observed that after certain value of <S2, the decrease of 
the sum of squares of difference is very small with the 
increase of Ci2. The point at which the decrease of the sum of 
the squares of the difference becomes very small is 
identified and the corresponding oZ value is taken as the 
accepted value. 


4.5.3 Calibration of 

When all the factors namely roadwidth, horizontal curvature 
and surface roughness are taken into consideration, the 
expression for Q reduces to 


In the above expression only 03 and Q values are 
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unknown. To calibrate 0!3 values for all the vehicle types we 

assume a positive value for 013 and the corresponding Q value 

is computed. Using this Q value V_. values for the cells 

3x 

FSC122 , FSC123 , FSC124 , FSS125, FSC132, FSC133 , FSC134 , FSS135 , 
FSC221 ,FSC222,FSC223,FSC224, FSS225,FSC231 , FSC232. FSC233, 

FSC234 , FSS235as shown in the Table 4.4.2 are predicted. The 
predicted values are then compared with the field values and 
their total sum of squares squares of the differences is 
computed. A graph is plotted between e<3 and the total sum of 
squares of the differences as shown in the Fig. 4. 5. 3. From 
the graph it is observed that at a particular value of C43, 
the decrease of the sum of the squares of difference becomes 
very small and that value is taken as 

The calibrated values of q , q , q , a , oi oi and 

X 2 >5 X 2 

Q values Ambassador car have been shown in the Table 4.2. 

This chapter explained in detail the various 
steps involved in calibration of Q value. The plot between 
field distribution and predicted distribution for Ambassador 
car in plain terrain with the calibrated Q values have been 
shown in Fig. 4.1 to 4.14. 

Next chapter explains the various experiments 
that have been carried out to validate the simulation model. 
The comparison between simulated and field data have been 


done in that chapter. 



TABLE 4.2 COMPUTED Q VALUES FOR DIFFERENT ROAD SECTIONS 


VEH.TYPE: AMBASSADOR TERRAIN: PLAIN 


ql= -1.96 

q2= -0.60 q3= 

-0.58 


1.0 

012= 2.0 0(3= 2 

.50 


SINGLE LANE ROADWAY 




Q112= -0.8443244 

0122= -0.7101527 

0132= 

-0.7303935 

0113= -0.9362803 

0123= -0.8111144 

0133= 

-0.7374662 

0114= -1.139228 

0124= -0.8443119 

0134= 

-0.8443119 

0115= -1.960000 

0125= -0.8906376 

0135= 

-0.7550169 

INTERMEDIATE LANE ROADWAY 

0211= -0.7572280 

0221= -0.7247202 

0231 = 

-0.6891277 

0212= -0.7667383 

0222= -0.7305164 

q232= 

-0.6861406 

0213= -0.8355168 

0223= -0.7514471 

0233= 

-0.6938245 

0214= -1.0015040 

0224= -0.7810684 

0234= 

-0.7033269 

0215= -1.9600000 

0225= -0.8266272 

0235= 

-0.7154139 


TWO LANE ROADWAY 


0311 = 

-0.6000000 

0321= 

-0.5954270 

0331 = 

-0.5915509 

0312= 

-0.6000000 

0322= 

-0.5949986 

0332= 

-0.5900677 

0313= 

-0.6000000 

0323= 

-0.5927517 

0333= 

-0.5882329 
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CHAPTER V 


VALIDATION OF THE SIMULATION MODEL 


5. 1 General 

Validation involves the verification that the 
structure of the model is correct and robust with respect to 
the parameter estimation. Generally the validation of a model 
is based on an agreement between the output of the model and 
available information from real life situation if they 
exist. The close such an agreement, better is the validity 
of the model. The model can be used to simulate future 
conditions only after it has been validated and found to 
replicate the dynamics of the process observed in the 
real world system under study. Validation thus 
provides an objective measure for confirmation of the 
correctness of the submodels along with assumptions involved 
in framing the model. Submodels like desired speed 
distributions, overtaking gap acceptance probability etc. 
that are based on the observed data must be validated 
separately. Using the overall simulation results, it is 
necessary to calibrate the values of some of the parameters 
and decision logics. The overall validation of the simulation 
model is a test of how well the submodels have been assembled 
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Validation involves the verification that the 
structure of the model is correct and robust with respect to 
the parameter estimation. Generally the validation of a model 
is based on an agreement between the output of the model and 
available information from real life situation if they 
exist- The close such an agreement, better is the validity 
of the model. The model can be used to simulate future 
conditions only after it has been validated and found to 
replicate the dynamics of the process observed in the 
real world system under study. Validation thus 
provides an objective measure for confirmation of the 
correctness of the submodels along with assumptions involved 
in framing the model. Submodels like desired speed 
distributions, overtaking gap acceptance probability etc. 
that are based on the observed data must be validated 
separately. Using the overall simulation results, it is 
necessary to calibrate the values of some of the parameters 
and decision logics. The overall validation of the simulation 
model is a test of how well the submodels have been assembled 



into a realistic representation of the system. 


of Ef f ootivorio&s f'or- Valiiiation 

The model outputs to be used as measure of 
ef f ectiver)ess in validation should be such that 

i) They represent the output of the whole system 
ii) They are capable of accurate measurements in the 
field. 

Significant disagreement implies that the model is 
unsatisfactory. We therefore decided to concentrate on 
analysis of the difference between the observed and simulated 
means and standard deviations in order to quantify just how 
"good" or "bad" the comparisons were. 

5.3 Test for Significance Of Difference Between Simulated 
aiul Observed Values 

Normal statistical procedure specifies the 
application of Chi square test for significance of difference 
between the two sets of data. This test is, however 
considered too stringent for model validation, as intuitively 
very good fits may also indicate ’significant difference’ 
when subjected to the test. Analysis of difference between 
observed and simulated values for mean values and standard 



“f r©6 spcsd Mas adopt.Bd as a cjuanti'tati.v© 
measure for the purpose. 

5. 4 Road Test Section for Model Validation ; 

For the purpose of validation, a hypothetical 
road stretch of length 4.5 Km and of road block length 200m 
each were generated. The different geometric feature that are 
being considered for validation are road width, horizontal 
curvature and surface roughness. Since the validation of the 
model have been carried out for free speed conditions, the 
traffic is generated only in one direction, so that movement 
of the vehicle is not constraint by the vehicles travelling 
in the other direction. 

5.4.1 Classification According To Roadwidth 

In road width category three different types of road 
width were considered. The three different classifications are 


Classification 

Name 

Road width 

I 

Single Lane 

3.5m 

II 

Intermediate Lane 

5.5m 

III 

Two Lane 

7.0m 



5.4.2 Classif icati'on According To Horizontal Curvature 


All the horizontal curves have been classified 
into four groups. The classification is done as follows. 

Classification Radius (R) 

I R <= 150m 

II 150m < R. <= 300m 

III 300m < R <= 450m 

IV 450m < R <= 1000m 

5.4*3 Classification According To Surface Roughness 

As far as surface roughness is concerned, the road 
has been classified into three different groups. First 

classification consists of all roads having surface 
roughness less or equal to 3000 mm/km. In the second 
classification, all the road having roughness between 3000 
and 6000 mm/km , are included. All the roads with 

roughness more than 6000 mm/km are in in the third group. 

For the each cell in the measurement matrix, a 
road file is created satisfying all the geometric 
combinations of that cell -The input given to the road 
submodel are, blockborder start coordinate, slope, horizontal 
curvature, roughness, roadwidth for each block. When road 



submodel creates a road file, it computes and assigns median 
block speed(V^^) and Q value for each vehicle, on each 
road block of the road stretch to be simulated. The Table 

5.1 shows a typical road file created by the road submodel. 
The Table 5.2 also shows the traffic file given as input to 
the simulation model. The power mass ratio used for different 
vehicles are listed in Table 5.3. Thus road files created for 

I 

each cells of the measurement matrix were given as input to 
the simulation programme along with the traffic file. Since 
execution time of the simulation program increases with the 
increase in vehicle number in the traffic file, we considered 
only seven vehicle types (Ambassador , Fiat, Maruti, Jeep, 
Bus.LCV and Truck), and .110 vehicles of each type were 
generated for each simulation run. After each simulation 
run, output file (event file) obtained was processed in the 
statistical submodel and results were tabulated. 

5.5 Comparison Of Simulated and Observed Data 

Comparison of simulated and observed data for single 
lane, intermediate lane and two lane roads for different 
horizontal curvature and roughness is presented in Table 

5.1.1 to Table 5.7.3. The analysis of data presented in 
the tables indicates good agreement between observed and 



vehicles studied. The percentage difference of the standard 
deviation of the mean free speeds Mere observed to be 
slightly high. The reason could be that, in the simulation 
model speed of a vehicle travelling on a particular road 
block remains constant. Whereas in the field, the speed does 
not remain constant throughout the block length. The model 
also gives slightly higher value of the of mean free speed 
as the radius of horizontal curvature decreases. This is due 
to the fact that Q value calibrated for cells in the 
measurement matrix having lesser horizontal curvature does 
not rotate these distributions as much as it should have 
beenCRef. Fig. 4.1). As a result free speed distributions 
predicted with these Q value are not very accurate. But in 
case of larger radius the rotation is perfect and predicted 
free speed distribution is very close to the field 
distribution. This is because the standard deviation of the 
mean free speed does not decrease as much as it should have 
been as mean free speed decreases. Thus this discrepancy of 
simulated data and field data on the road stretch having 
lesser radius of horizontal curvature can be attributed to 
errors field data collection. 

But as a whole, the percentage difference 
between the mean free speed for all the vehicle types and for 
all the road combinations were found to be less than 5%, 



parameters calibrated are correct. 


Thus it can be concluded that model developed is 
capable of predicting accurately the free speed distribution 
of vehicles on any road stretches. Next chapter presents the 
summary of the present Mork, conclusions and scope for the 
further work. 



TABLE 5.1 


ROAD FILE GIVEN AS INPUT TO THE SIMULATION PROGR/S 


Roadnidth = 7.0i, Radius of curvature = 300a, Roughness = 1500 la/ka 


4 0 1500 

4 700 1000 

4 1200 1180 

4 2044 900 

4 2900 1200 

4 3450 780 

4 4060 1300 

4 4510 900 

1 0 200 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 

0 30 0 0 0 1 

1 200 400 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 

0 30 1 1 0 1 

1 400 600 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 

0 14.35 -0.71 15.03 -0.87 17.21 -0.99 13.09 -0.64 12.50 -1.23 12.62 -0.60 14.61 -0.94 

0 11.76 -0.43 4.17 -0.42 2.22 -0.40 4.17 -0.40 8.34 -0.42 13.90 -0.41 9.73 -0.41 


0 30 1 1 0 1 



table 52 TRAFFIC FILE GIVEN AS INPUT TO THE SIMULATION PROGRAMME 


Identity 

No. 

Start 

Cord. 

Dest. 

Cord- 

Veh. 

Type 

Veh. 

BDS 

Cross 

BDS 

P/m 

Ratio 

Entry 

Time 

Entry 

Speed 

1 

0 

4500 

1 

18.47 

18.47 

10.00 

25.5 

15.7 

2 

0 

4500 

8 

14.85 

14.85 

10.18 

38.9 

12.6 

3 

0 

4500 

8 

14.62 

14.62 

2.26 

57.9 

12.4 

4 

0 

4500 

4 

16.68 

16.68 

7.94 

68.4 

14.2 

5 

0 

4500 

4 

16.10 

16.10 

2.50 

95.8 

13.7 

6 

0 

4500 

7 

19.81 

19.81 

3.70 

104.0 

16.8 

7 

0 

4500 

4 

16.09 

16.09 

5.51 

185.5 

13.7 

8 

0 

4500 

8 

18.89 

18.89 

2.64 

204.6 

16.1 

9 

0 

4500 

9 

15.35 

15.35 

3.46 

207.0 

13.0 

10 

0 

4500 

2 

18.19 

18.19 

4.76 

213.6 

15.5 

11 

0 

4500 

4 

17.38 

17.38 

12.70 

218.4 

14.8 

12 

0 

4500 

9 

16.91 

16.91 

4.27 

222.8 

14-4 

13 

0 

4500 

9 

16.56 

16.56 

2.70 

245.1 

14.1 

14 

0 

4500 

1 

18.47 

18.47 

10.00 

262.6 

15.7 

15 

0 

4500 

2 

17.02 

17.02 

8.40 

306,5 

14.5 

16 

0 

4500 

3 

22.59 

22.59 

5.77 

317.0 

19.2 

17 

0 

4500 

2 

19.12 

19.12 

7.00 

327.0 

16.3 

18 

0 

4500 

3 

21.52 

21.52 

6.99 

366.3 

18,3 

19 

0 

4500 

7 

15.81 

15.81 

2.22 

417.1 

13.4 

20 

0 

4500 

7 

15.91 

15.91 

3.20 

436.2 

13.5 

21 

0 

4500 

7 

18.04 

18.04 

3.38 

446.8 

15.3 



I ABLl 63 POWER MASS RATIO USED IN THE SIMULATION MODEL 


POWER CUMULATIVE PERCENTILE 

TO MASS 


RATIO 

AMBASSADOR 

FIAT 

MARUTI 

JEEP 

BUS 

LCV 

TRUCK 

1.0 

1.67 

1.60 

0.00 

1.67 

0.00 

1.67 

1.76 

2.0 

3.90 

3.90 

0.00 

3.90 

0.90 

3.90 

11.76 

3.0 

8.90 

8.89 

0.00 

8.89 

13.30 

8.89 

35.29 

4.0 

16.67 

16.67 

0.00 

16.70 

61.67 

16.70 

55.29 

5.0 

31.10 

31 .10 

0.00 

31.10 

86.70 

31.10 

69.41 

6.0 

44.40 

44.40 

1.60 

44.40 

93.30 

44.40 

80.00 

7.0 

63.30 

63.30 

3.90 

63.30 

98.30 

63.30 

88.23 

8.0 

74.40 

74.40 

16.70 

74.40 

98.30 

74.40 

94.11 

9.0 

83.30 

83.30 

31.10 

83.30 

100.00 

83.30 

98.82 

10.0 

88.90 

88.90 

44.40 

88.90 

100.00 

88.90 

100.00 

11.0 

94.40 

94.40 

63.30 

94.40 

100.00 

94.40 

100.00 

12.0 

95.60 

95.60 

74.40 

95.60 

100.00 

95.60 

100.00 

13.0 

98.89 

98.90 

83.30 

98.90 

100.00 

98.90 

100.00 

14.0 

100.00 

100.00 

88.90 

100.00 

100.00 

100.00 

100.00 

15.0 

100.00 

100.00 

94.40 

100.00 

100.00 

100.00 

100.00 

16.0 

100.00 

100.00 

95.60 

100.00 

100.00 

100.00 

100.00 

17.0 

100.00 

100-00 

98.90 

100.00 

100.00 

100.00 

100.00 

18.0 

100 . 00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

19.0 

100 . 00 

100.00 

100.00 

100.00 

100.00 

100-00 

100.00 

20.0 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 



TABLE 5.1.1 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: AMBASSADOR CAR LANE: SINGLE (3.5m) 


NUMBER OF VEHICLE SIMULATED: 110 


ROUGHNESS: 1500 mm/ km 


HORIZONTAL CURVE 

FIELD 

^0 

DATA 

a 

0 

SIMULATED 

^1 

t DATA 

ct 

1 

^0 

100 

^0 

*^1 ‘^O 

lOO 

^0 

R <= 300m 

46.63 

3.78 

48.00 

3.20 

+2.93 

-15.46 

300M < R <= 450m 

50.93 

3.78 

51.80 

2.80 

+ 1.70 

-26.03 

450 < R <= 1000m 

55.23 

3.89 

54.50 

2.80 

-1.32 

-28.04 

STRAIGHT 

59.53 

3.85 

58.00 

3.10 

-2.57 

-19.64 

ROUGHNESS :4500 mm/ km 

R <= 300m 

43.10 

3.10 

44.60 

2.70 

+3.48 

-13.53 

300m < R <= 450m 

45.98 

3.41 

46.80 

2.80 

+ 1.78 

-18.00 

450m < R <= 1000m 

48.87 

3.64 

48.30 

2.70 

-1.17 

-25.90 

STRAIGHT 

51.75 

4.19 

51.40 

2-90 

-0.67 

-30 . 80 


ROUGHNESS: 7500 mm/km 


R <= 300m 

37.10 

2.676 

37.7 

1.9 

+1.62 

-28.90 

300m < R <= 450m 

39.39 

3.207 

40.7 

2.7 

+3.33 

-15.80 

450m < R <= 1000m 

41.68 

2-948 

41.3 

2.3 

-0.91 

-21.90 

STRAIGHT 

43,97 

3.99 

43.6 

2.7 

-0.84 

-23.30 


jj U INDICATES MEAN, O . O INDICATES STANDARD DEVIATION 
^ 0 ’ 01 



TABLE 5.1.2 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: AMBASSADOR CAR LANE : INTERMEDIATE C5.5in) 


NUMBER OF 

ROUGHNESS: 1500 mm/km 

VEHICLE 

SIMULATED: 

110 



HORIZONTAL CURVE 

FIELD 

^0 

DATA 

*^0 

SIMULATED DATA 

^ ^1 

^ r^'o 
100 

^0 

°’l" *^0 
100 

R <= 150m 

46.99 

3.57 

47.80 

2.30 

+1.73 

-35.57 

150m < R <= 300m 

47.92 

3.64 

47.60 

2.30 

- 0.66 

-36.80 

300M < R <= 450m 

52.41 

4.38 

52.00 

3.70 

-0.78 

-13.61 

450 < R <= 1000m 

56.91 

4.40 

56.00 

3.00 

-1.59 

-32.58 

STRAIGHT 

61.41 

4.30 

58.90 

3.90 

-4.08 

-9.40 


ROUGHNESS: 4500 mm/km 




R 

<= 150m 

44 « 06 

3.45 

46.10 

2.20 

+4.63 

-36.23 

150m 

< 

R 

<= 300m 

44.96 

3.52 

45.20 

2.70 

+0.53 

-23.39 

300m 

< 

R 

<= 450m 

48.14 

4.07 

47.70 

3.20 

-0.91 

-21.53 

450m 

< 

R 

<= 1000m 

51-31 

4.12 

51.10 

2.70 

-0.41 

-34.59 




STRAIGHT 

54.48 

3.68 

53.90 

2.80 

-1.06 

-24.01 


ROUGHNESS :7500 mm/ktn 


R <= 150m 

38.85 

3.16 

38.50 

1 .70 

-0.90 

-46.20 

150m < R <= 300m 

38.58 

3.13 

38.90 

2.10 

+0.82 

-33.07 

300m < R <= 450m 

41.57 

2.68 

41.90 

1 .90 

+0.79 

-29.26 

450m < R <= 1000m 

44.56 

3.40 

44.90 

2.50 

+0.76 

-26.57 

STRAIGHT 

47.55 

3.50 

47.10 

2.30 

-0.95 

-34.43 

INDICATES 

MEAN SPEED IN KMPH, 

%'^1 

INDICATES 

STANDARD 

DEVIATION . 



TABLE 5.1.3 


COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: AMBASSADOR CAR LANE:TW0(7.0m) 
NUMBER OF VEHICLE SIMULATED: 110 

ROUGHNESS: 1500 mtn/km 


HORIZONTAL 

CURVE 

FIELD 

^0 

DATA 

*^0 

SIMULATED 

DATA 

^1 

^r^o 
100 

^0 

^1 ^0 

100 

a 

0 

R < = 

150m 

50.21 

3.71 

51.60 

2.70 

- +2.76 

-27.20 

150m < R <= 

300m 

51.61 

3.81 

51.30 

2.60 

-0.60 

-31.90 

300M < R <= 

450m 

55.35 

3.97 

54.60 

2.90 

-1.35 

-27.08 

450 < R <= 

1000m 

59.09 

4.21 

57.70 

3.10 

-2.35 

-26.43 

STRAIGHT 

62.82 

4.54 

59.30 

3.70 

-5.60 

-18.62 


ROUGHNESS: 4500 tntn/km 




R 

<= 150m 

46.66 

3.46 

47.70 

2.40 

+2.22 

-30.63 


150m 

< 

R 

<= 300m 

47.99 

3.55 

47.20 

2.40 

-1.64 

-32.54 


300m 

< 

R 

<= 450m 

50.88 

4.07 

50.40 

3.20 

-0.94 

-21.45 


450m 

< 

R 

<= 1000m 

53.77 

3.81 

53.30 

2.90 

-0.87 

-23.88 





STRAIGHT 

56.67 

4.20 

55.80 

3.00 

-1.53 

-28 . 62 



ROUGHNESS: 7500 mm/km 



R <= 150m 

41.45 

2.67 

40.80 

1.80 

-1.56 

-32.58 

150m < 

R <= 300m 

40.90 

2.63 

41.10 

1.80 

+0.50 

-31.55 

300m < 

R <= 450m 

44.11 

3.35 

44.70 

2.70 

+1.34 

-19.45 

450m < 

R <= 1000m 

47.31 

3.71 

47-10 

2.40 

-0.44 

-35.44 


STRAIGHT 

50.51 

3.66 

50.20 

2.60 

-0.61 

-29.13 


INDICATES MEAN SPEED 

IN KMPH, 


STANDARD 

DEVIATION 




TABLE 5.3.1 


COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: MARUTI CAR LANE: SIN6LE(3.5m) 
NUMBER OF VEHICLE SIMULATED: 110 
ROUGHNESS: 1500 mm/ km 


HORIZONTAL CURVE 

FIELD 

^0 

DATA 

a 

0 

SIMULATED 

^1 

DATA 

^0 

100 

^0 

*^1 ^0 
100 

^0 

R <= 300m 

54.65 

4.93 

55.50 

3.40 

+1.55 

-31.10 

300M < R <= 450m 

59.25 

4.93 

58.30 

3.40 

-1.60 

-31.00 

450 < R <= 1000m 

63.85 

4.72 

63.20 

3.80 

-1.02 

-19.57 

STRAIGHT 

68.45 

5.77 

68.10 

4.70 

-0.51 

-18.65 

ROUGHNESS: 4500 mm/km 


• 





R <= 300m 

49.41 

4.36 

51.30 

3.60 

+3.82 

-17.60 

300m < R <= 450m 

53.44 

4.87 

52.20 

3.60 

-2.32 

-26.16 

450m < R <= 1000m 

57.47 

4.69 

57.80 

3.20 

+0.57 

-31.88 

STRAIGHT 

61.51 

4.86 

61.90 

3.60 

+0.65 

-25.92 


ROUGHNESS: 7500 mm/ km 


R <= 300m 

41.84 

3.33 

41.20 

2.80 

-1.53 

-16.06 

300m < R <= 450m 

46.08 

4.04 

44.70 

3.00 

-2.99 

-25.79 

450m < R <= 1000m 

50.32 

4.52 

50.90 

3.70 

+ 1.15 

-18.17 

STRAIGHT 

54.56 

4.90 

55.60 

3.20 

+1.91 

-34.70 


jLi , /J INDICATES MEAN SPEED, INDICATES STANDARD DEVIATION 


TABLE 5.3.2 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: MARUTI CAR LANE: INTERMEDIATE(5 . Sm) 


ROUGHNESS: 

NUMBER OF VEHICLE 

15000 mm/ km 

SIMULATED 

: 110 



HORIZONTAL 

CURVE 

FIELD 

^0 

DATA SIMULATED 

O' u 

0 ^1 

DATA 

a 

1 

^0 

^XlOO 

^0 

‘^l ^0 
^XIOO 

% 

R < = 

150m 

56.48 

4.30 

56.90 

3.60 

+0.74 

-16.27 

150m < R <= 

300m 

57.20 

4.35 

56.90 

3.60 

-0.52 

-17.33 

300M < R <= 

450m 

61.80 

4.87 

62.20 

3.50 

+0.64 

-28.19 

450 < R <= 

1000m 

66.39 

5.58 

66.80 

4.10 

+0.62 

-26.56 

STRAIGHT 

70.98 

5.91 

70.50 

4.80 

-0.67 

-18.83 

ROUGHNESS: 

4500 mm/ km 






R < = 

150m 

52.59 

4.550 

55.0 

3.4 

+4.58 

-25.27 

150m < R <= 

300m 

53.37 

4.617 

53.8 

3.2 

+0.80 

-30.69 

300m < R <= 

450m 

56.97 

4.629 

57.4 

3.3 

+0.75 

-28.71 

450m < R <= 

1000m 

60.56 

4.704 

60.8 

3.6 

+0.40 

-23.46 

STRAIGHT 

64.16 

4.749 

64.5 

3.9 

+0.53 

-17.87 

ROUGHNESS 

:7500 mm/ km 




• 


R < = 

150m 

44.90 

3.940 

44.1 

2.3 

-1.12 

-41.62 

150m < R < 

= 300m 

44.17 

3.875 

44.8 

3.1 

+1.42 

-20.00 

300m < R < 

= 450m 

48.56 

1.133 

49.1 

2.6 

+1.11 

-37.09 

450m < R <= 

1000m 

52.95 

4.551 

53.2 

3.0 

+0.47 

-34.00 

STRAIGHT 

57.34 

4.400 

57.8 

3.3 

+0.80 

-25.00 


INDICATES MEAN SPEED, INDICATES STANDARD DEVIATION 




TABLE 5.3.3 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: MARUTI CAR LANE:TW0(7.0tn) 
NUMBER OF VEHICLE SIMULATED: 110 

ROUGHNESS: 1500 mm/ km 


HORIZONTAL CURVE 

FIELD DATA 

U c 

0 0 

SIMULATED 

^1 

DATA 

^0 

^XlOO 

% 

XlOO 

^0 

R <= 150m 

60.98 

4.66 

62.80 

4.50 

+2.98 

-3.43 

150m < R <= 300m 

61.92 

4.31 

62.30 

3.60 

+0.61 

-23.90 

300M < R <= 450m 

65.57 

5.11 

65.90 

3.60 

+0.50 

-21.72 

450 < R <= 1000m 

69.11 

5.74 

69.80 

4.40 

+0.83 

-23.43 

STRAIGHT 

72.87 

5.96 

72.20 

4.90 

-0.92 

-17.85 

ROUGHNESS: 4500 mm/km 

R <= 150m 

54.32 

4.88 

55.70 

3.50 

+2.54 

-28.27 

150m < R <= 300m 

55.92 

5.02 

56.20 

3.20 

+0.50 

-36.29 

300m < R <= 450m 

59.54 

4.54 

59.80 

3.50 

+0.43 

-22,92 

450m < R <= 1000m 

63.17 

4.75 

63.40 

3.80 

+0.36 

-20.03 

STRAIGHT 

66.79 

5.59 

66.90 

4.20 

+0.26 

-24.91 

ROUGHNESS: 7500 mm/ km 

R <= 150m 

48.41 

4.26 

47.90 

2.70 

-1.05 

-36.62 

150m < R <= 300m 

48.11 

4.23 

48.70 

3.20 

+1.22 

-24.40 

300m < R <= 450m 

52.31 

4.37 

52.70 

3.10 

+0.74 

-28.19 

450m < R <= 1000m 

56.51 

4.59 

56.80 

3.20 

+0.51 

-29.96 

STRAIGHT 

60.71 

5.01 

61.00 

3.60 

+0.50 

-26.00 


u . u INDICATES MEAN SPEED, 
^0 1 


, £>•, INDICATES STANDARD DEVIATION 
0 1 




TABLE 5.7.1 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: TRUCK LANE: SIN6LE(3.5m) 
NUMBER OF VEHICLE SIMULATED: 110 


ROUGHNESS: 1500 mm/km 


HORIZONTAL CURVE 

FIELD DATA 

0 0 

SIMULATED DATA 

^1 ^1 

/j - H a - a 

^ 1^0 10 
^XlOO ^XIOO 

^0 ^0 

R <= ’300m 

39.83 

3.26 

38.4 

2.80 

-3.59 

-14.21 

300M < R <= 450m 

45.75 

3.26 

45.0 

2.60 

-1.63 

-20.34 

450 < R <= 1000m 

51.66 

3.89 

50.3 

3.20 

-2.63 

-17.86 

STRAIGHT 

57.57 

4.73 

55.3 

4.10 

-3.94 

-13.37 


ROUGHNESS: 4500 mm/km 


R <= 300m 

38.72 

3.18 

39.20 

2.50 

+1.23 

-21.40 


300m < R <= 450m 

42.39 

3.07 

41.90 

2.40 

-1.15 

-22,05 


450m < R <= 1000m 

46.06 

3.19 

45.30 

2.70 

-1.65 

-15.36 


STRAIGHT 

49.74 

3.82 

48.60 

3.00 

-2.29 

-21.56 


ROUGHNESS: 7500 

mm/ km 







R <= 300m 

32.86 

2.64 

32.50 

1.90 

-1.95 

-28.24 


300m < R <= 450m 

35.87 

3.13 

35.40 

2.80 

-1.31 

-10.68 


450m < R <= 1000m 

38.89 

2.81 

38.40 

2.10 

-1.25 

-25.40 


STRAIGHT 

41.91 

3.01 

41.40 

2.40 

-1.21 

-20.42 



U , U INDICATES MEAN SPEED, O'., INDICATES STANDARD DEVIATION 

0 1 




TABLE 5.7.2 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: TRUCK LANE: INTERMEDIATECS.Stn) 
NUMBER OF VEHICLE SIMULATED: 110 
ROUGHNESS: 1500 mtn/km 


HORIZONTAL CURVE 

FIELD DATA 

0 0 

SIMULATED 

^1 

DATA 

^1 

^'r ^0 

^XlOO 

^0 

'^0 

^KlOO 

R <= 150m 

43.04 

3.24 

41.90 

2.40 

-2.65 

-25.92 

150m < R <= 300m 

42.84 

3.22 

42.00 

2.40 

-1.96 

-25.50 

300M < R <= 450m 

48.28 

4.05 

47.30 

2.80 

-2.03 

-30.91 

450 < R <= 1000m 

53.73 

4.45 

52.10 

3.60 

-3.03 

-19.13 

STRAIGHT 

59.18 

4.28 

56.50 

4.40 

-4.52 

+2.71 

ROUGHNESS: 4500 

mm/ km 






R <= 150m 

40.37 

3.36 

41.00 

2.30 

+ 1.56 

-31.54 

150m < R <= 300m 

40.66 

3.38 

40.20 

2.30 

-1.13 

-32.01 

300m < R <= 450m 

44.22 

3.65 

43.50 

2.50 

-1.62 

-31.60 

450m < R <= 1000m 

47.78 

3.76 

46.80 

2.80 

-2.05 

-25.60 

STRAIGHT 

51.34 

3.79 

50.00 

3.20 

-2.61 

-15.72 


ROUGHNESS: 7500 mm/ km 



R < = 

150m 

35.31 

3.07 

34.40 

1.80 

-2.58 

-41.36 


150m 

< R < = 

300m 

34.74 

3.02 

34.20 

2.20 

-1.55 

-27.15 


300m 

< R <= 

450m 

37.66 

2.81 

37.20 

2.10 

-1.22 

-25.34 


450m 

< R < = 

1000m 

40.58 

3.26 

40.10 

2.20 

-1.18 

-32.50 



STRAIGHT 

43.50 

3.32 

42.80 

2.50 

-1.61 

-24.78 



U , U INDICATES MEAN SPEED, O' , O' INDICATES STANDARD DEVIATION 
'^O 1 O 1 




TABLE 5.7.3 COMPARISON OF SIMULATED AND FIELD DATA 

VEHICLE TYPE: TRUCK LANE: TW0(7.0m) 
NUMBER OF VEHICLE SIMULATED: 110 
ROUGHNESS: 1500 mm/ km 


HORIZONTAL CURVE 

FIELD DATS 

*'0 % 

SIMULATED DATA 

^1 ^1 

^0 

^XlOO 

^0 

^1~ % 

^XIOO 


R <= 150m 

45.77 

3.17 

45.40 

2.70 

-0.81 

-14.80 


150m < R <= 300m 

46.24 

3.20 

45.50 

2.20 

-1.60 

-15.67 


300M < R <= 450m 

50.96 

3.57 

49.70 

3.20 

-2.47 

-19.13 


450 < R <= 1000m 

55.67 

3.99 

53.80 

3.90 

-3.35 

- 2.40 


STRAIGHT 

60.38 

4.16 

57.40 

4.50 

-4.93 

+8.12 


ROUGHNESS: 4500 mm/ km 

R <= 150m 

42.67 

3.10 

42.60 

2.50 

-0.35 

-17.64 


150m < R <= 300m 

43.36 

3.15 

42.90 

2.50 

-1.06 

-20.63 


300m < R <= 450m 

46.63 

3.57 

45.70 

2.70 

-1.99 

-24.55 


450m < R <= 1000m 

49.90 

3.81 

48.70 

3.00 

-2.40 

-21.30 


STRAIGHT 

53.17 

4.29 

51.70 

3.50 

-2.76 

-18.54 


ROUGHNESS: 7500 

mm/ km 







R <= 150m 

37.68 

2.79 

36.70 

2-00 

-2.60 

-28.31 


150m < R <= 300m 

36.79 

2.72 

37.00 

2.00 

+0.57 

-26.70 


300m < R <= 450m 

39.84 

3.20 

39.50 

2.20 

-0.85 

-31.33 


450m < R <= 1000m 

42.90 

2.82 

42.30 

2.50 

-1.39 

-11.37 


STRAIGHT 

45.96 

3.11 

45.20 

2.70 

-1.65 

-13.43 



Li Li INDICATES MEAN SPEED, O'., INDICATES STANDARD DEVIATION 
^0 * 1 ^ ^ 



CHAPTER VI 


SUMMARY CONCLUSION AND SCOPE FOR FURTHER WORK 

6. 1 Summary 

The primary objective of this thesis is to calibrate all 
the parameters involving free vehiclre model, form the field 
data collected by CRRI . The road submodel is then modified 
and finally traffic simulation model is extensively validated 
for free vehicle speed. 

As far as calibration of the parameters are 

concernrd, Q value which is the total transformation measure 

of basic median speed to the median speed V, , is 

Om 3m 

calibrated for ten different vehicle types. The total measure 

of rotation Q is also a function of q, ( ql for road width, 

q2 for horizontal curvature and q3 for surface roughness) 

with median speeds V, , V„ and V, included as weighing 

Im 2m 3m 

factor. Thus in order to calibrate Q value , q^^ values were 

calibrated first and then weighing factor a. values were 

calibrated. Using the calibrated Q value, distribution of the 
speed about the median value termed as speed 

distribution of the desired speed were generated. The 

generated distributions were then compared with actual 



distribution of the desired speed were generated. 


The 


generated distributions were then compared with actual 
distribution, and found satisfactory. 

Once the Q value is calibrated, it is implemented into 
the road submodel of the simulation model. Then hypothetical 
road stretch of 4.5 km. long, satisfying all the geometric 
combinations of each and every cell of the measurement matrix 
(as shown in Table 3.1) is generated. The simulation program 
is then run using the generated road as input and simultion 
result is then compared with field data. The percentage 
difference between simulated and field data were found to be 
less than 5% for each cell. 


6.2 Conclusion 

From the calibrated model we can evaluate the mean free 
speed of each vehicle type when moving under free condition. 

The results show a clear trend of decrease in the mean 
free speed as the road width decreases. The results also 
show that as the radius of horizontal curvature decreases, 
the mean speed decreases. The mean speed also decreases as 
the surface roughness increases. It is also observed that the 



vehicle type like Maruti, Mhich has higher power mass ratio 
has higher mean speed than the vehicle with lower power mass 
ratio like trucks. Hence it can be concluded that the results 
are consistent with respect to calibrated parameters. 

6.3 Scope for Further Work 

The traffic simulation model presented herein have 
covered much needed new ground in respect of both traffic and 
road parameters, and may be taken to mark an effective 
begining in simulation modeling for heterogeneous traffic on 
narrow road condition. Further work that . may be considered 
are 


1. Calibration and validation of interaction model 

2. Intersections, involving deceleration and 
acceleration of vehicles. 

3. Signalized intersections, involving the 
halt-proceed operations for vehicles - 


4. 


Bus-stop obstructions on narrow roads. 



can be 


5. Shoulder width, type and quality 
explici tlytaken into account and their effect on 

free speed can be obtained using the similar 
procedures given in this thesis. In this respect 
it is also worthwhile to predict the effect of 
wet shoulder as against dry shoulder. 

6. Environmantal conditions such as the inclemental 
weather conditions can be also studied. 

7. Free speeds were measured only during day time on 
the sample of roads used in this study. Since a 
large percentage of vehicles namely trucks and 
buses are found to operate at night it is 
desirable to study’ the free speeds of these 
vehicles at night. This is all the more important 
as there are a number of accedients taking place 
during night. There are numerous roadway factors 
affecting the speed of vehicles' at night and. for 
effecient traffic operations it is important to 
identify these factors and incorporate in design 



and maintenance phases of the roads. To achive 


this night vision cameras may be required 
inbuilt infrared optical system. 


Mith 
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